Carbon Nanotube with Square Cross-section: An Ab Initio Investigation 
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Recently, Lagos et al. (Nature Nanotechnology 4. 149 (2009)) reported the discovery of the 
smallest possible silver square cross-section nanotube. A natural question is whether similar carbon 
nanotubes can exist. In this work we report ab initio results for the structural, stability and 
electronic properties for such hypothetical structures. Our results show that stable (or at least 
metastable) structures are possible with metallic properties. They also show that these structures 
can be obtained by a direct interconversion from SWNT(2,2). Large finite cubane-like oligomers, 
topologically related to these new tubes were also investigated. 



I. INTRODUCTION 

The study of the mechanical properties of nanoscale 
systems presents new theoretical and experimental 
challenges^ 2 -. The arrangements of atoms at nano and 
macro scales can be quite different and affect electronic 
and mechanical properties. Of particular interest are the 
structures that do not exist at macroscale but can be 
formed (at least as metastable states) at nanoscale, spe- 
cially when significant stress/strain is present. Examples 
of these cases are atomic suspended chains 3 - - — and helical 
nanowires^. 

More recently^, it was discovered the smallest metal 
(silver) nanotube from high resolution transmission elec- 
tron microscopy (HRTEM) experiments. These tubes 
are spontaneously formed during the elongation of sil- 
ver nanocontacts. Ab initio theoretical modeling! sug- 
gested that the formation of these hollow structures re- 
quires a combination of minimum size and high gradient 
stress. This might explain why these structures had not 
been discovered before in spite of many years of investi- 
gations. Even from theoretical point of view, where low 
stress regimes and small structures have been the usual 
approach, no study predicted their existence. The unex- 
pected discovery of this new family of nanotubes suggests 
that such other 'exotic' nanostructures may exist. One 
question that naturally arises is whether carbon-based 
similar nanotubes (i.e., carbon nanotubes with square 
cross-section - CNTSC) could exist (figure [l}. 

From the topological point of view CNTSC tubes 
would require carbon atoms arranged in multiple square- 
like configurations. Molecular motifs satisfying these 
conditions, the so-called cubanes, do exist and they are 
stable at room temperature (figure [2]^. Cubane (CsHg) 
is a hydrocarbon molecule consisting of 8 carbon atoms 
in an almost perfect cubic arrangement (thus its name). 
Hydrogen atoms are bonded to each carbon atom (fig- 
ure [2]) , completing its four-bond valency. The 90 degrees 
angles among carbon atoms create a very high strained 
structure. Because of this unusual carbon hybridiza- 



tion, cubane was considered a 'platonic' hydrocarbon 
and believed to be almost impossible to be synthesized^. 
However, in 1964 Eaton and Cole^ achieved its break- 
through synthesis. Since then the cubane chemistry 
evolved quickly^kiS. Solid cubane^ proved to be remark- 
ably stable and polymers containing up to substituted 40 
cubanes units have been already synthesized^. However, 
up to now no tubular structure has been reporte d 12 ' 13 . 
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FIG. 1: (color online). Structural nanotube models. Frontal 
and lateral views: (a) Silver tubei, (b) SWNT(2,2), and (c) 
CNTSC. See text for discussions. 
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In this work we have theoretically investigated struc- 
tural, stability and electronic properties of CNTSC 
tubes. We have considered infinite (considering peri- 
odic boundary conditions) and finite (oligomers up to 
ten square units, figure [2]) structures. 

II. METHODOLOGY 

We have carried out ab initio total energy calculations 
in the framework of density functional theory (DFT) , as 
implemented in the DMol 3 code^. Exchange and corre- 
lation terms were treated within the generalized gradient 
(GGA) functional by Perdew, Burke, and Ernzerhofi^. 
Core electrons were treated in a non-relativistic all elec- 
tron implementation of the potential. A double numer- 
ical quality basis set with polarization function (DNP) 
was considered, with a real space cutoff of 3.7 A. The 
tolerances of energy, gradient, and displacement conver- 
gence were 0.00027 eV, 0.054 cV/A and 0.005 A, respec- 
tively. Initially we optimized the CNTSC unit cell with 



TABLE I: DMol 3 results for crystalline carbon allotropic 
structures: (a) Graphite, (b) Cubic diamond, (c) SWNT(2,2), 
and (d) CNTSC. 





(a) 


(b) 


(c) 


(d) 


Energy/atom (Ha) 


-38.085 


-38.081 


-38.041 


-38.003 


Lattice parameters: 










a (A) 


2.46 


3.57 


20.0 


20.0 


c (A) 


6.80 


3.57 


2.53 


1.62 


C-C bond-length (A): 


1.423 


1.537 


1.447 


1.580 








1.464 


1.616 



ultra-small carbon nanotubes (USCNTs) have been the- 
oretically investigated^— , it remains the smallest CNT 
experimentally observed with an estimated diameter of 
3 A 21 ' 22 . Graphite and diamond were also included in 
our study, because they are the two most stable carbon 
form and to provide a benchmark for the relative stability 
between the different tubes and these structures. 





FIG. 2: (color online). Cubane and its 'polymerized' units. 
The label refers to the number of square cross-sections in the 
structure. We considered structures from cubane up to n=10. 
Stick models, C and H atoms are in grey and white colors, 
respectively. 



fixed a and b parameters set to 20 A in order to ensure 
isolated (non-interacting) structures. The axial c lattice 
parameter was varied, and the total energy per unit cell 
calculated. All internal atomic coordinates were free to 
vary. Total energy versus unit cell volume was fitted us- 
ing the Murnaghan procedure^ to obtain the equilibrium 
c lattice parameter. For comparison purposes we have 
also considered graphite, diamond and carbon nanotubc 
SWNT(2,2). SWNT(2,2) was chosen because, although 



III. RESULTS AND DISCUSSIONS 

The results are presented in table |TJ As expected, 
graphite is the structure with the lowest energy (most 
stable), followed by diamond, SWNT(2,2) and CNTSC, 
respectively. Although the CNTSC energy per atom 
is relatively high (in part due to the strained C-C 
bonds, as in cubanes) its relative energy difference to 
SWNT(2,2) (0.0384 Ha) is similar to the difference be- 
tween SWNT(2,2) and cubic diamond (0.0395 Ha). 

In figure [3] we present the binding energy per unit 
cell. As a direct comparison it is not possible because 
the number of atoms in the minimum unit cell is differ- 
ent for SWNT(2,2) and CNTSC (eight and four, respec- 
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O 

Lattice parameter (A) 
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FIG. 3: Binding energy per unit cell as a function of ax- 
ial c lattice parameter for SWNT(2,2) (circles) and CNTSC 
(squares). It is also shown its interconversion curve (trian- 
gles). See text for discussions. 
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FIG. 4: (color online). Snapshots from the axial compres- 
sion process, showing the interconversion of SWNT(2,2) to 
CNTSC. (a) Initial SWNT(2,2), (b) and (c) intermediates, 
and (c) final CNTSC structures. 



there is a pathway that permits a direct interconversion 
from SWNT(2,2) to CNTSC. 

In figure [4] we present a sequence of snapshots from 
the simulations representing the interconversion process. 
The strain energy injected into the system by the radial 
compression (figure HJa)) produces a c-lattice expansion, 
leading to a structural transition (figure Hid)). The com- 
pression process produces new C-C bonds followed by 
carbon rehybridizations. The processes is better visual- 
ized in the supplementary materials (videol)^. 

We then proceeded with the CNTSC electronic struc- 
ture analysis. In figure [5] we displayed the band struc- 
ture and the density of states (DOS) for SWNT(2,2) and 
CNTSC tubes. Both structures present a finite DOS at 
the Fermi energy, characteristic of metallic regimes. Al- 
though CNTSC exhibits non-usual carbon hybridization 
it follows the general trends that small diameter carbon 
nanotubes are metalli c 24 ' 25 . 

It is possible that synthetic methods used to produce 
CNTs (such as laser ablation, chemical vapor deposition 
and arc discharge^) could also produce CNTSC, specially 
inside CNTs of different chiralitics, as in the case of 
SWNT(2,2)2i^. Another possibility could be a poly- 
meric synthetic approach, like the topochcmical ones to 
produce carbon nanotube of specific types that have been 
recently discussed in the literatur o 26 ' 27 . Considering that 
cubane molecules and their polymers exist and are sta- 
ble, we decided to investigate the local stability and en- 
dothermicity of cubane-like tubular oligomers that are 
topologically related to CNTSC. We carried out DMol 3 
calculations for the molecular structures displayed in fig- 



tively), we used a double CNTSC unit cell. As can be 
seen from figure El the results suggest that stable (or at 
least metastable, as a well defined minimum is present) 
CNTSC structures are possible. 

Our results also suggest that a direct interconversion 
from SWNT(2,2) to CNTSC is unlikely to occur via ax- 
ial(longitudinal) stretching. The extrapolation of the 
stretched SWNT(2,2) curve (figured! circle data points) 
could be misleading suggesting that it would be possible 
an interception with the stretched CNTSC curve (figure 
[3j square data points). However, this did not occur, the 
SWNT(2,2) can not preserve its tubular topology when 
its c- value is beyond 3.2 A. 

We then investigated whether if an assisted intercon- 
version would be possible, in our case we considered a 
continuously decrease of the tube radii value (in order 
to mimic an applied external (radially) pressure) while 
keeping the tube free to expand/contract longitudinally 
(see figure H] and video I— ). 

We have performed these calculations starting from an 
optimized SWNT(2,2) unit cell and then continuously 
decreasing its radii value and re-equilibrating the system 
and measuring the new c- values (figure El triangle data 
points). Our results show that under these conditions 



(a) SWNT(2,2) 




r Z 40 80 120 

Density of States 



(b) CNTSC 




r Z 20 40 60 80 

Density of States 



FIG. 5: Band structure and total density of states (in elec- 
trons/Ha) results for the (a) SWNT(2,2) and (b) carbon 
square-cross-section CNTSC. Energy is relative to Fermi level 
(dashed horizontal lines). Primitive unit cells have eight and 
four carbon atoms for SWNT(2,2) and CNTSC, respectively. 
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TABLE II: Total energy per carbon atom (in Ha) for the 
structures shown in figure [2] The corresponding value for the 
infinite structure is also presented. 



n 


e t (n) 


2 


-38.647 


3 


-38.429 




QC QOQ 
-OD.OiO 


5 


-38.259 


6 


-38.216 


7 


-38.186 


8 


-38.163 


9 


-38.145 


10 


-38.131 


CNTSC 


-38.003 



ure [2J The terminal C atoms were passivated with H 
atoms. 

In table [TT| we present the results for the total en- 
ergy per carbon atom as a function of the number of 
square cross-sections. Our results show that although the 
oligomer formation would require an cndothcrmal (ener- 
getically assisted) process, the structures are stable and 



the energy per carbon atom converges assyntotically to 
the corresponding value of the infinite tube (see also sup- 
plementary materials^) . 

In summary, based on a recent discovery of the small- 
est possible silver nanotube with a square cross-section 7 -, 
we have investigated whether a similar carbon-based 
structure could exist. We have used ab initio method- 
ology to investigate the structural, stability and elec- 
tronic properties of carbon nanotubes with square cross- 
section (CNTSC). Our results show that stable (or at 
least metastable) CNTSC (finite and infinite) structures 
can exist. They also show that it is possible to con- 
vert SWNT(2,2) to CNTSC under radial contraction. 
CNTSCs should share most of the general features ex- 
hibited by 'standard' nanotubes. Although the CNTSCs 
have not yet been observed, we believe our results had 
proven their feasibility. We hope the present work will es- 
timulate further works for these new family of carbon 
nanotubes. 
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